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Abstract. 

We demonstrate that characteristic nodules emerge in cultures of several malignant 
pleural mesothelioma (MPM) cell lines. Instead of excessive local cell proliferation, 
MPM nodules arise in culture by Myosin Il-driven cell contractility. Accordingly, 
the aggregation process can be prevented or reversed by suitable pharmacological 
inhibitors of Myosin II activity. A cell-resolved elasto-plastic model of the multicellular 
patterning process predicts that the morphology and size of the nodules as well as the 
speed of their formation is determined by the mechanical tension cells exert on their 
neighbors, and the stability of cell-substrate adhesion complexes. For small tension 
forces nodules are slow to develop and localized at the boundary of cell free areas. In 
contrast, a high intralayer tension quickly transforms all cell-covered areas into dense 
clusters interconnected by multicellular strands. Model simulations also indicate that 
a decreased stability of cell-substrate adhesions favors the formation of fewer, but 
larger clusters. Linear stability analysis of a homogenous, self-tensioned Maxwell fluid 
indicates the unconditional presence of the patterning instability. 
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1. Introduction 

Cell contractility and forces exerted on the cell’s microenvironment constitute an 
important mechanism of multicellular patterning. Morphogenetic movements in 
epithelia are particularly well known to utilize cell contractility: for example, an 
anisotropic contractile activity gives rise to cell intercalation, a process altering cell 
connectivity/adjacency in such a way that the whole tissue elongates in one direction 
while narrows along the perpendicular direction nia. Similarly, contractility-driven 
constriction of the free (apical) surface of the epithelium can give rise to bending 
or budding within an epithelial sheet | 3 ]. Mesenchymal cells can also contract the 
surrounding extracellular matrix (ECM): if a cell aggregate is placed on the surface of 
a collagen gel, cell traction reorganizes the collagen and create bundles of ECM that 
radiate away from the aggregate HE]. This observation led to the development of 
the mechano-chemical theory of pattern formation 012], according to which cells exert 
traction forces on an underlying deformable substrate and the resulting strain transports 
(convects) both cells and the ECM. Furthermore, strain-oriented ECM hlaments can 
guide cell motility, as cells are more likely to move parallel with the orientation of 
the ECM [HI El [in] . This mechanism was suggested to guide vascular patterning, and 
endothelial cells were reported to be able to detect and respond to substrate strains 
created by the traction stresses of neighboring cells m- 

Cellular contractility is also an important factor in tumor progression as it can 
contribute both to local and distant spreading of malignant cells [HI OSl IH [ig. Here 
we focus on the role of cell contractility in malignant pleural mesothelioma (MPM), a 
tumor arising from mesothelial cells lining the pleural cavity. This highly aggressive 
disease has an incidence of 1 in 100 000 in Europe [I6]. Despite recent advances in the 
treatment of MPM, it has an extremely poor prognosis with almost all patients dying 
from their tumor. A characteristic, pathognomonic feature of MPM is the formation of 
multiple, macroscopic pleural tumor nodules, which - due to the special two dimensional 
environment - may pinch off and contribute to the local spreading of malignant cells. 

In this manuscript we demonstrate that several MPM cell lines can form nodule- 
like aggregates in vitro when cultured at high cell density. Time-lapse analysis of their 
formation and experiments with myosin Il-specihc inhibitors reveal that cell contractility 
is a key factor in the nodule formation process. To interpret these experimental 
hndings, we propose a cell-resolved elasto-plastic model of a contractile cell layer. By 
computational simulations we explore the effects of key model parameters, such as the 
the stability of cell-cell and cell-substrate adhesion and the magnitude of the cell-exerted 
contractile forces. 
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2. Methods 

2.1. Cell lines 

SPClll cells were established from human biphasic MPMs and kindly provided by 
Prof. R. Stahel (University of Zurich, Switzerland). P31 cells were a kind gift from 
Prof. K. Grankvist (University of Umea, Sweden). The Meso cell lines were established 
by the Vienna MPM group as described recently im. NP3 normal mesothelial cells 
were isolated from pneumothorax patients using the same protocol [T8] . 

2.2. Culture conditions 

Gelatin-coated dishes were obtained by incubating 10% gelatin solution-B (Sigma) in 
PBS for 45 min at room temperature. We also prepared 1.7 mg/ml Gohagen-I (Gorning) 
gels according to the manufacturer’s instructions. 

Gells were grown at 37°G in a humidihed, 5% GO 2 , 95% air atmosphere. 
The DMEM (Lonza) medium was supplemented with 10% FGS (Invitrogen) and 
1% penicihin-streptomycin-amphotericin B (Lonza). MPM cells were cultured on 
plastic tissue culture substrates (Greiner). NP3 cells are, however, more sensitive to 
environmental factors and have a limited proliferative potential. Thus, NP3 cells were 
cultured in gelatin-coated dishes, for up to three passages. 

2.3. Reagents 

To interfere with normal Myosin II function, we utilized Y27632, the rho kinase 
(ROGK) inhibitor (Merck Millipore) and Blebbistatin (Merck Millipore), an inhibitor of 
actomyosin crosslinking. Y27632 was solved in destilled water, stored in the form of a 10 
mM stock solution, and used at 100 fiM hnal concentrations in DMEM. Blebbistatin was 
solved in DMSO, stored as 50 mM aliquots, and used at 40-80 fiM hnal concentrations. 

2 . 4 . Immunostaining 

SPGlll cells were seeded at conhuent density (40000 cells/cm^) either on the surface 
of collagen gels or on glass coverslips and were grown for 6 days. Samples were hxed 
using 4% PEA for 15 min at 4‘’G, permeabilized with 0.25% Triton-XlOO for 10 min 
at 4°G, and incubated with the following antibodies: Polyclonal anti-hbronectin (1:100, 
Ab2033, Millipore); polyclonal anti-beta-catenin (1:100, G2206, Sigma). After washing, 
sections were incubated for 30 min with anti-rabbit Alexa-488 secondary antibody 
(Life Technologies, Garlsbad, GA). Filamentous actin was stained by Phalloidin-TRITG 
(P1951, Sigma). Samples were analyzed by confocal laser scanning microscopy using 
the Bio-Rad MRG-1024 system (Bio-Rad, Richmond, GA). 
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2.5. Physical sections 

Physical cross sections were obtained from cultnres grown on the snrface of collagen-I 
gels. Samples were embedded in Spnrr’s mixtnre. Semithin (1 fim) sections were cnt 
perpendicularly to the surface of the collagen and stained with 0.5% Toluidin blue. 
Images were taken using Zeiss Axioskop 2 microscope equipped with a lOOx objective 
and coupled to an Olympus DP50 camera. 


2.6. Time-lapse microscopy 

Time-lapse recordings were performed on a computer-controlled Leica DM IRB inverted 
microscope equipped with a Marzhauser SCAN-IM powered stage and a lOx N-PLAN 
objective with 0.25 numerical aperture and 17.6 mm working distance. The microscope 
was coupled to an Olympus DP70 color CCD camera. Cell cultures were kept at 37C 
in humidihed 5% CO 2 atmosphere during imaging. Phase contrast images of cells were 
collected consecutively every 10 minutes from each of the microscopic helds. 

To maintain high cell density cultures for several days, we 3D printed three 6 mm 
diameter mini-wells into regular 35 mm tissue culture dishes. The side wall of the wells 
was formed by fused polylactic acid hlaments. This setup allowed us to achive high 
(10®/cm^) cell density with a medium/cell ratio characteristic for culture with ten times 
lower cell density (10^/cm^). The 3 ml medium was replaced every 5 days. Four modihed 
culture dishes, each containing 3 wells, were observed by time-lapse videomicroscopy for 
up to 14 days. We recorded four adjacent helds from each well, thus we collected 48 
images in each imaging cycle (every 10 minutes). 


2.1. Computational model 


To explore the dynamics of cell contraction-driven aggregation and nodule formation, 
we adapted our cell-resolved elasto-plastic computational framework [19] into a two 
dimensional model. Briehy, forces distributed along the membranes of mechanically 
coupled adherent cells are replaced by a single net force and a torque, simplifying the 
system to a network of particles and elastic beams (Fig. [^. 

In the computational model torques and shear forces arise due to relative movement 
of adjacent particles. Hook’s law determines the force needed to uniaxially compress 
or stretch cells (Fig. Id). In particular, Fj^ is associated with changing the length of 
link I which connects particles i and j located at r* and rj, respectively: 


Fi =k{\rj-ri\-ii). ( 1 ) 

In Eq. ^ ii is the equilibrium length of link I and fc > 0 is a model parameter, 
characterizing the stiffness of the cytoskelon. 

A link I exerts a torque if its preferred direction at particle i, ti^i, is distinct from 
its actual end-to-end direction (Figs. [^,d). We assume that the torque exerted on 
particle i is proportional to the difference between the preferred and actual directions: 


Mi,/ = X Uij), 


( 2 ) 
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Figure 1. Mechanical model of multicellular clusters. Links represent 
mechanically connected cytoskeletal structures of adherent cells. Rectangular 
shapes indicate cell membranes, not explicitly resolved in the model, a: 
Two particles, i and j interconnected with a link in a mechanical stress-free 
configuration. Unit vectors, tj and t^, which co-rotate with the particles, 
represent the preferred direction of the link, b: Compressed cells. The 
interaction of the two cells gives rise to spatially distributed forces (red arrows), 
which are replaced by the repulsive net forces Fj and Fj. c: A symmetric 
rotation of both particles yields torques Mj and Mj acting on particles i and j, 
respectively. These torque vectors are perpendicular to the plane of the figure, 
d: A lateral misalignment of the particles creates torques M*, Mj and also 
shear forces Fj, Fj acting at the particles. 


where the microscopic bending rigidity g, a model parameter, can be calibrated from 
macroscopic material properties of the tissue [12]. The condition for mechanical 
equilibrium, together with Eqs ([^ and ([^ allow the calculation of the forces and torques 
within the system [T9] . 

Multicellular plasticity is modeled using rules that rearrange the network of 
intercellular adhesions. The stability of intercellular adhesion complexes depends on 
the tensile forces they transmit [20] • Thus, in our model the probability of removing 
link I during a short time interval At follows Bell’s rule [21] as 

PiAt = At, (3) 

where is a threshold value, and 

f ; = e(b')d' (4) 

is the tensile component of the force transmitted by the link and A is a scaling factor 
which characterizes the stability of connections. 

Mechanical connections can be established between two Voronoi neighbor particles, 
i and j. We assume, that during a short time interval At the probability of inserting a 
new link is a decreasing function of the distance dij between the particles: 

qijAt = B (l - At. (5) 

\ ^max J 
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The scaling factor B represents the intensity of cellnlar protrnsive activity devoted to 
scanning the environment and the ability to form intercellnlar contacts. The maximal 
distance cells explore for new connections is denoted by dmax- 

Simnlations are event-driven: nsing the probability distribntions (|^ and ([^, we 
generate the next event /i and waiting time r according to the stochastic Gillespie 
algorithm izg. The waiting time nntil the next event is chosen from the distribntion 

logP(r) =-r + J , (6) 

V I i,j ) 

where the sums are evaluated by iterating over existing links I as well as over all possible 
Voronoi neighbor particle pairs i,j not connected by a link. 

2.8. Simulation Parameters 

Two dimensional initial conditions were generated by randomly positioning N = 400 
particles in a square of size L = 20. The unit distance of the simulations was set to 
the average cell size, do ~10/im, thus the 2D cell density is 1 cell/unit area. In the 
initial condition we enforced that the distance of two adjacent particles is greater than 
dmin = 0.8do. Particles that are Voronoi neighbors are connected by links when their 
distance is less than dmax = 2do. For a mechanical stress-free initial configuration we 
set the ti^i preferred link direction vectors as well as the equilibrium link lengths ii so 
that no internal forces or torques are exerted in the system. As a boundary condition, 
we fixed the position and orientation of the particles located near the perimeter of the 
simulation domain. 

The mean waiting time between simulation events is set by parameters A and B. 
We choose our time unit as 1/B ^ 10 min, the time needed for two adjacent cells 
to establish a mechanical link. We set the lifetime of an unloaded link to 1/B 1 

day. Thus, according to these values, two cells pulled away by a force separate in 
~ 10 h, a characteristic value consistent with the time scale observed in our cell culture 
experiments. 

3. Results 

3.1. MPM cells form nodules in vitro 

To study the long-term behavior of mesothelioma monolayers, we cultured several, 
human patient-derived MPM cell lines for up to a few weeks. In this time frame MPM 
cells form macroscopic multicellular aggregates or nodules, which can be 50 pm thick and 
reach a millimeter in lateral extent (Fig. [^. The time needed to generate nodules varies 
across the cell lines. When SPClll and p31 cells are seeded at confluency, aggregates 
reach a macroscopic size after 3-5 days. In contrast, nodule formation in the Meso53, 
Meso62 and MesoSO lines takes a few weeks. Formation of similar nodules is rather 
uncommon in cultures of other epithelial cells. Most importantly, under similar culture 
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Figure 2. Mesothelioma cells spontaneously form nodules in vitro. Panels a-e 
depict five distinct human MPM cell lines, SPClll, Meso62, p31, MesoSO and 
Meso53, after one (a, c, e) or two (b, d) weeks in culture. In comparison, 
NP3 human primary mesothelial cells remain in a monolayer and do not form 
aggregates even after 17 days in vitro (f). Asterisks mark nodules, arrow 
indicates strands interconnecting nodules. 


conditions and dnration NP3 human primary mesothelial cells remain in a monolayer 
(Fig-i). Long-term NP3 cultures reach a quiescent state in which both the number of 
cell divisions and cell deaths are substantially reduced. 

Physical cross-sections of the nodules reveal densely packed cells containing 
prominent, pleiomorphic nuclei. Immunohistochemistry with antibodies against the 
extracellular matrix protein hbronectin indicates the presence, but not the abundance 
of ECM in the MPM nodules. Thus, the in vitro MPM nodules are dense, highly cellular 
structures. 



Contractile tumor nodules 



Figure 3. Time course of spontaneous nodule formation in a culture of SPClll 
cells. Panels a-e are frames from a time-lapse recording taken at seeding, 24h, 
48ti, 60h and 70h after seeding, respectively. The initial confluent monolayer 
(a) develops cell density fluctuations (b, c), which eventually become three 
dimensional nodules interconnected by multicellular strands (d-e). f: The 
aggregation process is visualized as a kymogram, where pixels along the same 
line are plotted for each frame. Lines from earlier frames are located at the top 
of the image. Scale bars correspond to 150 /rm. 
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Figure 4. In vitro MPM nodules are rich in stress filaments and are 
mechanically integrated. a: Actin filaments are visualized by confocal 
microscopy using TRITC-phalloidin (red). Nodules are interconnected with 
multicellular strands that exhibit parallel bundles of stress filaments, b: 
Stress cables within multicellular strands align across cell membranes (arrows, 
visualized by catenin antibodies - green). 


Time-lapse microscopy recordings of the aggregate formation process (Fig. reveal 
that nodules do not form by unusually active cell proliferation. Instead, cells move 
towards nodules or interconnecting strands. Kymograms (pixels along a selected line 
plotted for several consecutive frames) indicate a gradual decrease of the extent of an 
aggregate, a decrease which is approximately linear in in time. 

3.2. Nodule formation is driven by actomyosin contractility 

Visualization of actin filaments with fluorescent phalloidin (Fig.|^ reveals the abundance 
of stress cables in the strands, organized into parallel bundles connecting the adjacent 
nodules. Higher resolution confocal images indicate that several stress hlaments reach 
across cell bodies and even form structures that are continuous across cell membranes. 

Both the observed cell movements and the presence of profound stress cables within 
the aggregates suggest that Myosin 11 dependent cell contractility is an important driving 
force to collect MPM cells into nodules. To test this hypothesis, we administered 
drugs that interfere with normal Myosin 11 activity. Blebbistatin is a potent blocker 
of acto-myosin contractility as it reduces the affinity of myosin heads to actin [23] . The 
compound Y27632 is a specific inhibitor of Rho-associated kinase (ROCK), one of the 
activators of Myosin. 

Both blebbistatin (data not shown) and Y27632 (Fig. could substantially hinder 
or completely eliminate both the stress cables and nodule formation in cultures of 
SPClll cells. When inhibitors were administered at high concentrations (100 /iM 
for Y27632 and 80 /iM for blebbistatin), the previously formed nodules flattened: 
their height decreased by up to 50% while simultaneously extended laterally. When 
the inhibitor is present from the onset of the culture, nodule formation can be fully 
prevented (data not shown). The inhibitors are reversible: aggregation resumes after 
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3DIV 


3 DIV, 

5 h treatment 


7 DIV, 

4 days treatment 


10 DIV, 

3 days after wash 


Figure 5. In vitro nodule formation requires Myosin II activity. Frames of a 
time-lapse recording show the morphology of two parallel cultures, an untreated 
control (a,c,e,g) and one which was transiently exposed to Y27632, a Rho- 
associated kinase inhibitor (b,d,f,h). The two cultures are similar at 3 days 
in vitro, at the onset of Y27632 treatment (a,b). Five hours long exposure 
to Y27632 is sufficient to induce simultaneous flattening and lateral extension 
of the nodules (d). Cell-dense three dimensional structures, clearly visible in 
untreated control cultures (e), are largely absent after 4 days of ROCK inhibitor 
treatment (f). Removal of the drug restores the contractility of MPM cells: 
three days after replacing the medium nodule morphologies are similar in both 
cultures (g,h). 
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replacing the medium with fresh DMEM, and three days after medium replacement 
nodule morphologies are similar in the unperturbed and previously myosin-blocked 
cultures. 


3.3. Computational model of contractile cells 


To obtain a computational model of contractile cell sheets, we augmented our elasto- 
plastic particle model with two new rules: one that controls equilibrium link lengths to 
maintain a steady intercellular tension, and another to represent adhesions between the 
cells (particles) and the culture substrate. 

We assume that cells strive for a specific contractile environment which they 
maintain as a homeostatic state: cells increase their contractility if tensile forces between 
adjacent cells are below a target value, F*. In our model this is achieved by reducing 
the equilibrium link length di in Eq. ([^ . Conversely, when tensile forces are too strong, 
the equilibrium length of the links is increased. In particular, we assume that the rate 
of change in the link length is proportional to the difference between F* and F/, the 
tensile component (|^ of the force transmitted by the link: 

ddi ^ , Ff - F* , , 

— Y* —(^) 


where 1/C ~ 1 h sets the temporal scale of the feedback regulation. The last term in ([^ 
is an uncorrelated (white) noise, representing random cell shape changes due to factors 
not considered explicitly in our model. Equation was integrated over a time interval 
At, elapsed until the next stochastic event effecting particle connectivity, resulting: 


di(t + At) = dAt) -\- Cdo—!— - At + a/ DAtf, (8) 

where ^ is a pseudo-random variable with unit standard deviation and parameter 
D ^ 1/xm^/h controls the noise amplitude. 

The mechanical equilibrium of surface-attached cells requires the net force exerted 
by adjacent cells to be balanced by the force transmitted through the cell-substrate 
adhesion complexes to the substrate. Thus, in our model for each particle i we introduce 
the equilibrium position r° and angle 0°. The force (and torque) transmitted to the 
substrate is given by the linear relations 


Fi,o = fco(r* - r°) (9) 

and 

Mi,o = go{(j)i - 0 -)e 2 . ( 10 ) 


As we discuss below, external forces acting on a cell may contribute to its displacement 
trough complex and mostly unexplored processes. Here we assume that cells tend to 
move in the direction of the net external force, hence the eqilibrium position is updated 
according to 


dr? ado 
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Figure 6. Time development of a typical simulation, at t=2.5 h (a), t=7.5 h 
(b), t=13 h (c) and t=16 h (d). Blue to green colors indicate increasing tensile 
stress within the links. Simulation parameters are listed in table I. 


and 


d^ 

dt 


/S 




M 


i,0- 


( 12 ) 


Parameters a and (3 set the external force-related bias in cell movements, while the noise 
in Eq. (1^ gives rise to random walk-like movements. 


3.4- Simulations 

Simnlations of the contractile monolayer (Fig. reveal a pattern formation seqnence 
which is in several aspects analogous to the process seen in vitro (Fig.|^. First, the slight 
initial inhomogeneity of the monolayer is amplihed resulting in the formation of holes. 
Then, cell free areas continue to grow as their boundary is unstable: links constituting 
the boundary need to balance the pulling forces exerted by the bulk of the cells. Thus, 
connections at boundaries experience higher tensile stress which increases both their 
length and the probability of their removal. Both effects extend the area of the holes. 
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t [h] t [h] 

Figure 7. The progress of patterning, characterized by the standard deviation 
of particle density, S, as a function of time. Pattern formation is faster for 
higher target force F* values (a) and for higher values of a, the motion bias 
towards external forces (b). 


As a consequence, cell density increases in the rest of the system, and eventually cell- 
covered areas appear as contractile nodes connected by linear strands. While nodules 
develop in several MPM cultures (Fig. |^, the appearance of cell-free areas is observable 
only in a subset of MPM lines. For example, SPClll or p31 nodules form on top of 
a basal monolayer and in such cases our model corresponds to the dynamic upper cell 
layers. 

During the aggregation process the spatial distribution of cells becomes more 
inhomogeneous. Thus, the progress of patterning can be characterized by S, the 
standard deviation of the coarse-grained particle density held. While the time course 
shown in Fig.j^is characteristic for all simulations, two model parameters have important 
roles in determining both the morphology and the speed of the aggregation (Fig. [^. One 
such parameter is the ratio between the steady state contractile force of the cells, F*, 
and the Bell threshold of adhesion stability, F^. Figure [^compares conhgurations that 
are at the same, late stage {S ~ 4) of the patterning process. For small forces {F* < F^) 
high cell density clusters develop at the boundary of cell free areas, similar to the in 
vitro patterns observed for the MesoSO and Meso53 lines (Fig. |^. In this case the 
aggregation process is slower - again, in accord with empirical in vitro data. For large 
forces {F* > F^) links with high particle density interconnect similarly dense nodules, 
reminiscent to the patterns observed in SPClll cultures. 

Model parameter a, characterizing the magnitude of external force-directed cell 
displacements, sets the spatial scale of the pattern (Fig. [^. When cell-substrate 
adhesions are stable (a <C 1), smaller clusters develop which are close to each other. In 
contrast, when cells respond strongly to external forces (i.e., substrate adhesion is weak 
or highly adaptable), fewer and larger clusters form. 
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Figure 8. Morphologies characteristic for various target force values. 
Simulations performed with F*/= A (a), F*/Fq = 2 (b), F*/F^ = 1 (c), and 
p* jpO = Q _5 a,re shown at the same stage of pattern formation (S' = 3.85). 
For small forces F* < F^ high cell density clusters develop at the boundary 
of cell free areas (arrowheads), while the particle density remains low far from 
such boundaries (asterisk). In contrast, for large forces F* > F^ the particle 
density is more uniformly high. 


3.5. Linear stability analysis 

To better understand the patterning mechanism we performed a linear stability analysis 
of the spatially homogenous state. As a continuum model of the multicellular system, 
we consider a viscoelastic Maxwell material which relaxes shear stresses through an 
exponential decay, similar to the behavior of the particle model [19]. In one dimension 
the stress distribution a{x, t) is related to the velocity held of the tissue, v{x, t), through 
the equation 


( 13 ) 
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Figure 9. Morphologies characteristic for various adhesion parameter values. 
Simulations performed with a = 0.3 (a), a = 0.1 (b), a = 0.03 (c), and a = 0.01 
(d) are shown at the same stage of pattern formation (S = 3.5). For smaller 
values of a, the size of the aggregates and the characteristic distance between 
aggregates (arrows) decreases. 


where E is the macroscopic elastic modulus, determined by model parameters k and g 
[19], and 1/C is the rate of macroscopic stress adjustment analogous to the parameter 
C in Eq. ([^. The hrst term on the right hand side represents elastic stress and the 
second term describes the relaxation of the active contractile stress generated by the 
cells to a value cT*(p) that increases with the local cell density as cr*(p) ~ F*p. 

Velocity is set by stress divergence as 

- ( 14 ) 


dx 


= OLV 


where the drag coefficient a is analogous to the a pararmeter in Eq. lED' 

In the absence of cell death and proliferation the cell density p(a;, f) satisfies the 
conservation equation 

dp d ^ r. 
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Considering no-flnx bonndary conditions the total cell nnmbers are conserved so that the 
average cell density po stays constant in time. Thns the spatially nniform eqnilibrinm 
solntion of the above system is p{x, t) = po, t) = 0, a{x, t) = ao = cr*(po). 

We can analyse the linear stability of the nniform state by considering the time 
evolntion of small pertnrbations of the form 

p(a;,t) =po + pe^*+*''" (16) 

v{x,t) (17) 

a{x,t) = + (18) 

where s is the growth rate of a monochromatic pertnrbation of wavennmber q, in the 
linearised system assnming that the pertnrbation amplitndes p,v,a are small. 

After snbstitntion and neglecting higher then linear terms we obtain 

sp + iqpov = 0 (19) 

sa =iEqv - C{a - a'^{po)p) (20) 

iqa = aa (21) 

Eliminating the amplitndes we obtain the following qnadratic eqnation for the growth 
rate s 


+ s 


C + Eq^- -q- 
a 


,<(P«)EC 


= 0 


a 


with the solntions 

s±(«) = ^ (-%) ± V%F-i-Va) 


( 22 ) 


(23) 


where 


a = 


<{po)EC 


h{q) = C + Eq^- 


(24) 

a a 

Since a and b{q) are positive, one of the solntions is always negative while the other 
is always positive for all wavennmbers. Thns the spatially nniform eqnilibrinm state is 
nnconditionally nnstable resnlting in increasing cell density flnctnations. This leads to 
the formation of cell aggregates, however as the pertnrbations grow the development of 
the system is not described by the linear approximation since the neglected nonlinearities 
become important. 


4. Discussion 

4.1. Cell displacements and external forces 

As in the case of cell-cell connections, mechanical load acting on cell-snbstrate adhesion 
complexes rednces their lifetime [20]. To relate cell movements and external forces acting 
on a cell we envision the following process: when a cell-snbstrate connection breaks, the 
same mechanical load is distribnted along the remaining adhesion complexes. Thns, 
each of the remaining adhesion sites transmits a larger force, their strain is increased 
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leading to a small displacement of the cell body in the direction of the net external force 
acting on the cell. Furthermore, when new adhesion complexes form, their equilibrium 
(stress-free) conhguration will correspond to the actual, slightly shifted position of the 
cell. Thus, by detaching and re-attaching adhesion complexes, the cell relaxes the shear 
stress between its cytoskeleton and the adhesion substrate, and moves in the direction 
of the external force. In addition to this purely mechanical connection, external stress 
may also effect the polarity of active cell migration [2H 1^ . 

4-2. Contraction and aggregation 

Aggregation and sorting involves acto-myosin contractility within the cortex of zebrahsh 
cells [2H]. Here we show that cell groups may also contract through stress cables spanning 
across several cells, and the resulting system self-organizes into expanding cell-free 
areas and eventually into free-standing aggregates. Similar behavior also takes place 
at much smaller length scales in the cytoskeleton of individual cells, where acto-myosin 
contractility gives rise to f-actin bundles. For example, a similar contractile system, 
but one that also includes diffusion, has been studied by [27] in the context of pattern 
formation on the actomyosin cell cortex in which contractility is regulated at molecular 
level. Further extensions of this problem have been described recently |2Sl 122] . 

4-3. Future treatment options 

The demonstrated ability of myosin-II inhibitors to flatten mesothelioma nodules may 
open a new therapeutic method. As cells within avascular, three dimensional nodules 
are less exposed to systemic drugs, efficient inhibition of nodule formation could enhance 
the effective drug concentration at the targeted cells. In this case, we would expect a 
synergistic effect between existing anti-cancer drugs and myosin-II inhibitors. 
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